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ABSTRACT 

Motivated by the discovery of prolate rotation of stars in Andromeda II, a dwarf spheroidal com¬ 
panion of M31, we study its origin via mergers of disky dwarf galaxies. We simulate merger events 
between two identical dwarfs changing the initial inclination of their disks with respect to the orbit 
and the amount of orbital angular momentum. On radial orbits the amount of prolate rotation in 
the merger remnants correlates strongly with the inclination of the disks and is well understood as 
due to the conservation of the angular momentum component of the disks along the merger axis. For 
non-radial orbits prolate rotation may still be produced if the orbital angular momentum is initially 
not much larger than the intrinsic angular momentum of the disks. The orbital structure of the 
remnants with significant rotation is dominated by box orbits in the center and long-axis tubes in 
the outer parts. The frequency analysis of stellar orbits in the plane perpendicular to the major axis 
reveals the presence of two families roughly corresponding to inner and outer long-axis tubes. The 
fraction of inner tubes is largest in the remnant forming from disks oriented most vertically initially 
and is responsible for the boxy shape of the galaxy. We conclude that prolate rotation results from 
mergers with a variety of initial conditions and no fine tuning is necessary to reproduce this feature. 
We compare the properties of our merger remnants to those of dwarfs resulting from the tidal stirring 
scenario and the data for Andromeda II. 

Subject headings: galaxies: dwarf - galaxies: Local Group - galaxies: fundamental parameters - 
galaxies: kinematics and dynamics - galaxies: structure 


I. INTRODUCTION 

The most widely recognized theory of structure for¬ 
mation in the Universe is based on the assumption of 
the existence of cold dark matter, significantly domi¬ 
nating the baryonic one. In such a Universe, galax¬ 
ies form by hierarchical merging of smaller systems. 
Large galaxies show plenty of signs of past or ongo¬ 
ing mergers but it seems harder to collect observational 
evidence on the mass scale of dwarf galaxies. Never¬ 
theless, some serious effort has already been made in 
order to reveal signs of galactic mergers down to ^ 
10^ (iGeha et al.ll2QQ5l: iMartmez-Delgado et al.l[2QI2l: 

iCrnoievic et al.ll2QI4HPandel et al.ll2QI5[). 

In parallel, attempts have been made to understand 
dwarf mergers from the point of view of numerical simu¬ 
lations. In this case high resolution is needed in order to 
adequately model a low-mass galactic merger in a cosmo¬ 
logical context because of the differences in scales. Al¬ 
though limited in resolution, studies in the context of the 
Local Group demonstrated that a significant fraction of 
present-day dwarfs underwent a sig nificant merger with 
another dwarf galaxy in the past (jKlimentowski et al.l 
l2QIQl : [Deason et M]l2QI4f ). These mergers most likely oc¬ 
cur early on, at the outskirts of the Lo cal Group and be¬ 
tween bound pairs or within groups. iKazantzidis et al.l 
(|2QII[ ) demonstrated that mergers between initially disky 
dwarf galaxies may produce pressure-supported systems 
with parameters very similar to those of dwarf spheroidal 
(dSph) galaxies found in the Local Group. These sys- 
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terns seem particularly suitable for testing the merger 
paradigm on small scales as their proximity offers a pos¬ 
sibility for high-resolution measurements. 

A few candidates for merger remnants have been al¬ 
ready identified among the dSph galaxies of the Lo¬ 
cal Group. The first was the Fornax dwarf for which 
shells of stars in the vicinity of the main body were 
observed (|Goleman et al.l I2QQ5D . Other cases may in¬ 
clude the most distant dSphs such as Getus and Tu- 
cana that could not have interacted with the Milky Way 
and MSI very strongly in the past and be tidally trans- 
for med so the merger scenar io seems particularly plausi¬ 
ble (|Kazantzidis et aTll2QIl[ ). though hypothetical as no 
direct observational evidence exists in this case. 

Recently, a new merger candidate appeared in the form 
of the Andromeda II (And II) dSph galaxy, one of the 
brightest satellite s of MSI. The o bject has a number of 
peculiar features. iHo et al.l (|2QI2[) detected strong rota¬ 
tion signal along the 2D minor axis which means that 
the object probably has a streaming motion around the 
longest axis in SD, i.e. prolate rotation. Note that And II 
is th e only dSph with a c l ear pr olate rotation reported so 
far. iMcGonnachie et~al] (|2QQ7[) identified a few distinct 
stellar populations in And II d iffering in age, m etallicity 
and density profiles. However, iHo et al.l (|2QI2f ) found no 
significant difference in the kinematics of the metal rich 
and metal poor stars. The presen ce of the metal licity 
gradient was recent ly confirmed by iHo et all (|2QI5f ) and 
I Vargas et al.l (|2QI4f ). I Vargas et al.l (|2QI4D also found an 
intermediate age population which makes And II an ex¬ 
ception among their sample of six fainter dSphs that pos¬ 
sess only populations of sta rs older than 10 Gyr. More¬ 
over, [Amorisco^^F^I] (|2Q14D reported the detection of re- 
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TABLE 1 

Initial conditions for the simulations. 


Merger 

Vx 

[kms-1] 

Vy 

[kms-i] 

Vz 

[kms-1] 

Ax/Atot 

LYlLtot 

LziLtot 

Line type 

I90R 

-8.0 

0.0 

0.0 

0.71 

0.0 

0.71 

red solid 


8.0 

0.0 

0.0 

0.71 

0.0 

-0.71 

red dashed 

I60R 

-8.0 

0.0 

0.0 

0.87 

0.0 

0.50 

cyan solid 


8.0 

0.0 

0.0 

0.87 

0.0 

-0.50 

cyan dashed 

I120R 

-8.0 

0.0 

0.0 

0.50 

0.0 

0.87 

magenta solid 


8.0 

0.0 

0.0 

0.50 

0.0 

-0.87 

magenta dashed 

I90VY1 

-8.0 

1.0 

0.0 

0.71 

0.0 

0.71 

green solid 


8.0 

-1.0 

0.0 

0.71 

0.0 

-0.71 

green dashed 

I90VY2 

-8.0 

2.0 

0.0 

0.71 

0.0 

0.71 

blue solid 


8.0 

-2.0 

0.0 

0.71 

0.0 

-0.71 

blue dashed 


Note. — Odd rows list the initial conditions for Dwarf 1 and even rows for Dwarf 2 of a given merger. The first column identifies the 
merger simulation (see text for details). The following columns list the initial components of orbital velocities and angular momenta (in 
units of the total angular momentum of each dwarf). The last column indicates the type and color of the line used for a given dwarf and 
simulation in the figures throughout the paper. 


g:i Qns with lower v elocity dispersion in the kinematic data 
of lHo et ^ (|2Q12[ ) which can be interpreted as a presence 
o f a stellar stream in And 11. 

iLokas et al.l (|2Q14a[ ) proposed a merger scenario for the 
origin of And 11. In this model, the merger occurs be¬ 
tween two dwarf galaxies initially composed of a disk 
and a dark matter halo. The dwarfs had equal masses 
and only differed by the disk scale lengths. They were 
placed on a radial orbit towards each other and their 
disks were inclined by ±45 deg to the collision axis so 
that the angular momentum vectors of the disks and the 
velocities of the dwarfs were initially in the same plane. 
The model adequately reproduces the observed proper¬ 
ties of And II, including the prolate rotation and shape. 
The presence of the different stellar populations can be 
interpreted as originating from different dwarfs: the two 
populations then have different density distributions (as 
a result of different initial disk scale lengths) but very 
similar kinematics. 

Another scenario known to lead to the formation of 
dSph galaxies is the one based on tidal interaction of the 
initially disky dwarf sat ellites with bigger hosts, like the 
Milky Way or M31 fe.g. iLokas et al.ll2Qi4bl . |2()15[) . How¬ 
ever, as argued bv iLokas et al.l (|2Q14a[) . tidal interaction 
is able to remove very efficiently the initial rotation of the 
disk (around the shortest axis) as it transforms to a bar 
and then a spheroid, but cannot induce any significant 
prolate rotation (around the longest axis). 

The initial condition s used in the simulation described 
in iLokas et al.l (|2Q14aD may seem rather special making 
the merger highly improbable. In this paper we aim to 
investigate how much these conditions may be relaxed 
and still reproduce the prolate rotation in the final prod¬ 
uct of the merger. We study the effect of different initial 
inclination angles of the dwarf disks and the effect of 
non-radial orbits of the dwarfs. In section 2 we present 
the simulations used in this study. The orbital evolution 
of the dwarfs is described in section 3 while their shapes 
and kinematics in section 4. Section 5 is devoted to a 
more detailed study of the orbital structure of the rem¬ 
nants. In section 6 we compare the predictions of the 
scenario based on mergers with an alternative way to 
form dSph galaxies that relies on their interaction with 
a bigger host. In section 7 we compare our predictions 


to the real data for And II and remark on the possible 
future, more complete model that will describe all the 
available data. The summary and concluding remarks 
follow in section 8. 

2. SIMULATIONS 

We performed five collisionless simulations of merg¬ 
ers between two identical dwarf galaxies. The simula¬ 
tions differed in the initial inclinations of the disks in the 
dwarfs and their relative velocities. The numerical real¬ 
ization of the dwarf galaxy is the sa me as used pr e viously 
in s imulations of tida l stirring by ILokas et al.l (j2Q14bf ) 
and lLokas et al.l (j2Q15l ). The model was generated f ollow - 
ing t he method des c ribed in IWidrow fc Dnbinskil (j2QQ5l ) 
and IWidrow et al.l (j2QQ8[ ) which allows to create self- 
consistent, equilibrium Wbody models o f galaxies with 
an exponential disk and an NEW halo (iNavarro et al.l 

[T^. 

Our dwarf galaxy had a dark matter component in 
the form of the NEW halo of total mass Mh = 10^ Mq 
and concentration parameter c = 20. The exponen¬ 
tial disk had the mass of Md = 2 x 10^ Mq, the scale- 
length i?d = 0.41 kpc and thickness z^/Rd = 0.2. Each 
component of the dwarf contained 10^ particles (which 
makes the total of 4 x 10^ particles per simulation). The 
mergers were followed using the Wbody simulation code 
GADGET-2 (|Springell 120051 ) with the adopted softening 
scale lengths of 0.02 kpc and 0.06 kpc for the disk and 
halo particles, respectively. Each merger was evolved 
for 10 Gyr with positions and velocities of the particles 
stored every 0.05 Gyr. 

We chose the coordinate system of the simulation box 
so that its origin was at the center of gravity of the 
two dwarfs. Dwarf 1 was always placed at coordinates 
(X, T, Z) = (—25, 0,0) kpc and Dwarf 2 at (25, 0,0) kpc. 
Eor radial mergers X is the original collision axis, for the 
other ones XY is the original collision plane. In order to 
investigate the effects of different disk inclinations and 
merger orbits, we keep the same model of mass distri¬ 
bution for both dwarfs in all simulations. We also keep 
the same initial distances, 50 kpc, and the same relative 
velocity in the direction of the X axis, 16kms“^. At 
the beginning of each simulation, the angular momentum 
vectors of the dwarf disks lie in the XZ plane. Of course. 
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Fig. 1.— Relative distances (top panel) and velocities (bottom 
panel) of the dwarfs for simulations with different initial relative 
velocities. 

due to the finite precision of the computation, the dwarfs 
can slightly deviate from a given axis or plane during the 
simulations. 

The components of the initial angular momenta of the 
dwarf disks and the relative initial velocities of the dwarfs 
for all five simulations are listed in Table m The labels 
assigned to the simulations (the first column of the Ta¬ 
ble) refer to the starting configuration: the number after 
‘T stands for the inclination angle (in degrees) between 
the angular momentum vectors of the dwarfs, ‘R’ means 
radial merger, and for non-radial mergers ‘VY’ indicates 
the velocity (in kms“^) along the Y axis. 

We use capital letters to denote relative velocities or 
distances of the two dwarfs in the simulation as well as 
the original coordinate system we defined for the simu¬ 
lation (X,Y,Z). Small letters {x^y^z) refer to the coor¬ 
dinate system associated with the principal axes of the 
components representing the luminous matter of either 
individual dwarfs (Dwarf 1 or Dwarf 2) or of the whole 
galaxy resulting from the merger (i.e. disk components 
of both dwarfs combined). In the latter system the ma¬ 
jor axis corresponds to the x axis, the intermediate one 
to y and the minor one to z. Similarly, other quantities 
related to this system are marked with small letters, e.g. 
Vx is the rotation around the major axis of the respective 
dwarf or galaxy, and such quantities always refer to just 
the luminous matter. 

For comparisons of the results we divide the simula¬ 
tions into two groups: (1) radial merger simulations with 
different initial inclinations of the disks (I90R, I60R and 
I120R) and (2) simulations with the same initial incli¬ 
nation but with different initial relative velocities (I90R, 
190VY1 and I90VY2). The simulation I90R is included in 
both groups for reference and is treated as default. This 
simulat ion is closest to the one discussed bv iLokas et al.l 
(|2Q14a[ ) and shown to reproduce well the properties of 
And II, except for the fact the we now use identical dwarf 
progenitors. 

3. ORBITS OF MERGING DWARFS 



X [kpc] 

Fig. 2.— Trajectories of dwarfs in simulations with different ini¬ 
tial relative velocities. 

Fig. [1] shows the evolution of relative distances and ve¬ 
locities of the dwarfs for simulations with different initial 
relative velocities. For these simulations, the trajectories 
of the dwarfs in the collision plane (the XY plane) are 
displayed in Fig. O For our radial mergers (along the X 
axis), the orbits are almost the same in all three cases. 

In the case of radial mergers the initial conditions and 
the subsequent evolution are very symmetrical. As a re¬ 
sult, the major axis of the merger remnant is aligned 
with the merger axis (the X axis of the simulation box). 
By introducing some velocity in the Y direction (for sim¬ 
ulations 190VYl and 190VY2) we introduce asymmetry 
in our simulations because the merger becomes prograde 
for Dwarf 1 and retrograde for Dwarf 2. This causes dif¬ 
ferences in the properties of the two dwarfs at the end 
of the simulation for non-radial mergers. All dwarfs ini¬ 
tially have the velocity of 8kms“^ towards the common 
center of mass. Simulations 190VYl and 190VY2 have 
additional 1 km s“^ and 2 km s“^, respectively, in the per¬ 
pendicular (Y) direction. The second approach of the 
dwarfs is on a close-to-radial orbit in both cases but for 
190VY2 it takes about 0.7 Gyr longer to reach it. 

For I90VY1 the merging times and results are quite 
similar to simulation I90R, although there is less rota¬ 
tion in the final product. The major axis of the merger 
remnant is in the XY plane but inclined by about 22 deg 
with respect to the X axis. For I90VY2 the major axis 
of the remnant also stays in the XY plane, now inclined 
by 67 deg to the X axis, but otherwise the simulation 
outcome is very different. There is just a little rotation 
in the final product and there are significant differences 
in the shape and rotation between Dwarf 1 and Dwarf 2 
at a given time step. For this reason, in the follow¬ 
ing figures showing the properties of the dwarfs, we plot 
just the results for Dwarf 1 except for simulation 190VY2 
for which we show results for Dwarf 2 as well (the blue 
dashed curves). 

Fig.[3]shows several snapshots from each of the simula¬ 
tions with different initial relative velocities. The dwarfs 
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Fig. 3. — Snapshots from the simulations with different initial relative velocities. The three rows of panels show dwarfs from simulation 
I90R, I90VY1 and I90VY2 respectively, in projection onto the orbital plane XY. Panels in different columns correspond to different stages 
of evolution, as marked in the upper left corner of each plot. Black lines indicate the orbits followed by the dwarfs up to a given time with 
the solid line for Dwarf 1 and the dashed one for Dwarf 2, respectively. The colors code the mean line-of-sight velocity. Each panel shows 
the area of 20 x 20 kpc. 


are shown in projection onto the orbital plane XY and 
the colors code the mean line-of-sight velocity. The first 
panel of each row shows the situation around the first 
pericenter passage. We note that Dwarf 1 of I90VY2 (on 
the right side of the panel) is notably more deformed than 
Dwarf 2. In the second-column panels dwarfs are located 
near their first apocenters and the third-column snap¬ 
shots roughly correspond to the time when the centers 
of the dwarfs have just merged, but the merger remnant 
is still far from equilibrium. The last-column panels are 
the final snapshots from the simulations, after 10 Gyr of 
evolution. 

4. KINEMATICS AND SHAPE 
4.1. Prolate rotation 

In order to check how much prolate rotation is present 
in each merger remnant, we found the principal axes of 
the stellar component of each dwarf separately and com¬ 
puted the mean rotation velocity around the major axis 
in the sphere of radius 2 kpc centered on the center of the 
respective dwarf. We plot the evolution of this rotation 
in time for Dwarf 1 for all simulations in Eig. 01 
Due to the symmetry of the initial conditions in radial 


mergers, the evolution of this rotation in Dwarf 1 is the 
same as for Dwarf 2 in the same simulation. Eor I90VY1, 
the asymmetry is too small to lead to significant differ¬ 
ences between the two dwarfs. On the other hand, in the 
case of I90VY2 the evolution is dramatically different for 
the two dwarfs, so we include the result for Dwarf 2 as a 
dashed curve in Eig. 01 

Eor all simulations, there is essentially no rotation 
around the intermediate axis and the evolution of the 
rotation around the minor axis is trivial: the disk pre¬ 
serves the rotation until the two dwarfs actually start to 
merge and then the rotation vanishes and remains close 
to zero until the end of the simulation. We show one 
example of the evolution of the mean rotation velocity 
around all principal axes for our default simulation I90R 
in Eig. [5l 

Eig. ini shows the maps of line-of-sight properties of the 
final outputs for all five simulations in projection onto 
the xy plane, i.e. as seen along the ^ axis of the merger 
remnant. Here, x, y and ^ correspond to the major, in¬ 
termediate and minor axis of the combined stellar com¬ 
ponents of both dwarfs. The images confirm that the ro¬ 
tation signal (middle-column panels of Eig. [6]) is around 
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time [Gyr] 

Fig. 4. — Evolution of the mean rotation velocity around the 
major axis for Dwarf 1 in the simulations with different initial rela¬ 
tive velocities (top panel) and for the radial mergers with different 
inclination of the disks (bottom panel). The blue dashed curve 
corresponds to Dwarf 2 from simulation 190VY2. 



Fig. 5. — Evolution of the mean rotation velocities around the 
principal axes of Dwarf 1 in simulation I90R. 

the longest (x) axis, i.e. the rotation is clearly prolate in 
all simulations, even when the rotation level is very low 
as in 190VY2. 

The amount of prolate rotation in the final product of 
a given simulation seems to be controlled in a simple way. 
For radial mergers, decreasing the initial angle between 
the angular momenta of the dwarfs and the merger axis 
(which corresponds to decreasing the angle between the 
angular momenta of the two dwarfs, from 120 to 90 and 
60 deg) we get more rotation in the final result. When 
the initial angular momenta of the disks are more aligned 
with the merger axis the components of the momenta 
along the merger axis are larger. Since these components 
are preserved during such symmetrical mergers the final 
products have more prolate rotation. 

For non-radial mergers, the more orbital velocity we 
add in the direction perpendicular to the radial, the less 
intrinsic rotation we get. This result can be interpreted 
by comparing the orbital and angular momenta of the 
dwarfs. In the case of I90VY1 (I90VY2) the orbital 
angular momentum of each dwarf is about twice (four 
times) larger than the intrinsic angular momentum of 
the stellar component of the dwarf. We conclude that 
the prolate rotation in the remnant cannot be produced 




log(i:) V|os [km/s] 0,05 [km/s] 


Fig. 6. — Maps of the surface density (first column), the mean 
line-of-sight velocity (middle column) and the line-of-sight velocity 
dispersion (last column) for the stars in the final outputs of the 
simulations as viewed along the minor axis of the combined stellar 
component. Different rows correspond to different simulations, as 
marked on the right. More vertical orientations of the initial disks 
and more radial merger orbits lead to stronger prolate rotation. 
The surface density of the stars is in units of 645 Mq kpc“^ with 
contours spaced by 0.5 in logS. The contours in velocity and 
velocity dispersion are spaced by 1 km s“^. 

if the orbital angular momentum during the merger is 
significantly larger than the intrinsic one. 

The maps of line-of-sight velocity dispersion in the 
right column of Fig.[6]show considerable variation among 
the different simulations without any clear trend: some 
remnants have a maximum dispersion in the center while 
others have the maxima at larger distances. It seems that 
the complicated structure of the dispersion maps may be 
used as an indicator of the merger origin of the galaxy 
even after many Gyr of evolution. 

4.2. Triaxiality 

In order to study the evolution of the shape of the 
merger remnants, we again analyzed each dwarf sepa¬ 
rately selecting the stars inside the sphere of 2 kpc cen- 
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Fig. 7. — Evolution of the axis ratios and the triaxiality param¬ 
eter for Dwarf 1 in different simulations. The gray horizontal lines 
in the lower two panels separate the regions of oblate, triaxial and 
prolate shape. The blue dashed curve corresponds to Dwarf 2 of 
simulation 190 VY2. 


I60R I120R 

I90R - I90VY1 



Fig. 8. — Position angles of the intermediate axes of the stellar 
components for the four simulations which exhibit figure rotation. 
The angle is measured in the YZ plane during the last 5 Gyr of 
each simulation. 

tered on a given dwarf. In Eig. [7] we plot the interme¬ 
diate to major axis ratio, b/a, the minor to major axis 
ratio, c/a, and the triaxiality parameter, T, defined as 
T = [1 — {b/a)^]/[l — (c/a)^]. When T is lower than 1/3 
the galaxy is rather oblate, values greater than 2/3 in¬ 
dicate a prolate shape and when T falls between these 
values the galaxy has a triaxial shape. 

The most dramatic change of the shape of the dwarfs 
happens around the second passage (approximately at 
^4Gyr). The ratio 6/a, which is equal to unity at the 
beginning, slightly drops to values between about 0.8 and 
0.9 for all simulations. The evolution of c/a is more in¬ 
teresting. It grows from the original value of 0.12 up to 
the value ranging from 0.43 for I90VY2 (Dwarf 1) to 0.86 
for I60R. Eor radial mergers, there is a simple rule for the 
behavior of c/a: the final ratio decreases when the ini¬ 
tial angle between the angular momenta of the dwarfs 
increases, i.e. the merger remnants become thinner. Eor 
merger 190VY2 b/a remains similar for both dwarfs but 
c/a is significantly different (0.43 for Dwarf 1 and 0.73 
for Dwarf 2). In the later stages of the simulations, all 
dwarfs can be classified as triaxial. Only I90VY1 crosses 
the line of the oblate shape in several time steps and 
Dwarf 2 of 190VY2 has more prolate values of triaxiality. 

Referring again to Eig. [6] we note that the shapes of the 
merger remnants also differ in terms of their boxiness. In 
this respect, the dependence seems to be the strongest 
for the radial mergers. Clearly, the simulation with the 
dwarf disks oriented most vertically (I60R) with respect 
to the merger axis the isodensity contours are most boxy. 
In order to make this statement more quantitative we cal¬ 
culated the Eourier mode ^44 of the surface distributions 
of the stars (such as those shown in the left column of 
Eig. [6]) for all simulations. The amplitude of the mth 
Eourier mode is given as = (1/Y) 
where (pj are the particle phases in the projection plane 
and N is the total number of particles. It turns out that 
over the last few Gyr of evolution A 4 ^ 0.06 for I60R 
while it is at least a factor of two lower for the other 
simulations. This shape must originate from the orbital 
structure of the merger remnant and we discuss this issue 
further in the following sections. 

5. ORBITAL STRUCTURE 
5.1. Figure rotation 
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Fig. 9.— Surface density of stellar particles in simulation I90R 
at two time steps: at the first encounter (upper row) and after the 
merger (lower row). Columns show projections along the Z, Y and 
X axis of the simulation box (from the left to the right). 


Although the ratios of principal axes do not vary sig¬ 
nificantly over the last 5 Gyr for all simulations, all but 
I90VY2 show clear figure rotation around the major x 
axis. The orientation of the major axis remains stable 
for the whole period of 5 Gyr and it coincides with the 
collision (X) axis for radial mergers and is inclined by 
22 deg from the collision axis for I90VY1 but remains in 
the XY plane. The evolution of the position angles of 
the intermediate axes for these four simulations are plot¬ 
ted in Fig. [HI For the purpose of these measurements 
we computed the principal axes for the combined stellar 
components from both dwarfs in each simulation. 

The periods of rotation of the intermediate (and short) 
axes are 1.7, 2.2, 6.0, and 3.7 Gyr for I60R, I90R, I120R, 
and 190VYl, respectively. The speed of the figure ro¬ 
tation is correlated with the rotation of stellar particles 
around the long axis (see Fig. SI) but the particles rotate 
at a higher speed (at least for the inner 2kpc). Also the 
direction of the particle rotation is the same as the figure 
rotation. 

The origin of the figure rotation can be traced to the 
tidal distortion of the disk just before the first encounter. 
We illustrate this in Fig. [9] which displays stellar density 
in simulation I90R in three projections along the Z, Y 
and X axis of the simulation box. The upper row shows 
the projections just before the first passage at 1.95 Gyr 
and the bottom row illustrates the situation soon after 
the two disks merged, at 5.3 Gyr. The upper left panel 
of the Figure (view along the Z axis) shows the distor¬ 
tion of the disks due to the tidal force acting before the 
merger. Since the stellar particles in both disks are mov¬ 
ing upwards in this plot in the parts of the disks that 
are already in touch, the tidal force distorted the disk so 
that the upper parts are closer than the lower ones. This 
results in the tumbling (clockwise around the X axis) of 
the elongated shape after the merger seen in the lower 
right panel (view along the X axis). 

Simulation I90VY2 has the major axis inclined by 
67 deg with respect to the X axis and lies in the XY 
plane also being stable for the last 5 Gyr of evolution. 
The position angles of the other axes vary by only about 
30 deg over the last 3 Gyr so this merger remnant does 
not show any significant figure rotation. 
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Fig. 10.— Fraction of a given orbit type as a function of radius for 
the stellar component in the final output of simulations I60R, I90R, 
I120R and I90VY1. Different line types correspond to different 
orbit types as shown in the legend. 

5.2. Orbit classification 

In this subsection, we focus on simulations which ex¬ 
hibit significant prolate rotation, so we exclude 190VY2 
from further analysis. 

At the beginning of the simulations stellar particles are 
situated in inclined disks on nearly circular orbits. Dur¬ 
ing the merger particle orbits are transformed to form 
an object with a shape that can be well approximated 
by a triaxial ellipsoid. In static triaxial potentials most 
particle orbits belong to one of a few major families: box 
orbits, short-axis tube orbits, an d inner and oute r long - 
axis tube orbits (|de Zeeuwlll985[) . iSchwarzschildl (|1982 ) 
also constructed a model of the triaxial ellipsoid with 
figure rotation around the short axis by properly pop¬ 
ulating different orbit families. Glassifying orbits in a 
non-static potential is more difficult. 

We performed a simplified classification of particle or¬ 
bits using the last 4 Gyr of the simulations when the 
galaxies are already merged. We defined one orienta¬ 
tion of the coordinate system {x^y'^z') which coincides 
with the actual principal axes of the galaxy only at one 
or few particular moments. For the rest of the time only 
the major axis is along the x coordinate. 

We ado pted modified criter ia of orbit classifica- 
tion from ISchwarzschildl (|1993D as recently used by 
iGaida et al ] (l2nT^ . Particles that keep a constant sign 
of the x-component of the angular momentum, are 
marked as long-axis tubes. Since the direction of the 
short axis is not stable, short-axis tubes are not classi¬ 
fied in our case. 

We also find orbits with low average values of the x, 
y' ^ and ^'-component of the angular momentum. Specif¬ 
ically, we require the sum of these average values, Lxy'z' i 
to be less than 0.1 Mq kpc^ Myr“^ (for comparison, in all 
outputs almost all particles inside 2 kpc have current to¬ 
tal angular momentum smaller than 1 Mq kpc^ Myr“^). 
When the figure rotation is pre sent the box orbits ac quire 
some net angular momentum (jSchwarzschildin~982[ ) thus 
we call this group just low angular momenta orbits. It 
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is worth noting that Ly/ and Lz' of most of the particles 
average out to zero and Lxy'z' is typically dominated by 
the average value of 

The rest of the particles, which do not fall into one of 
the above categories, are marked as other. These parti¬ 
cles have Lxy'z' > 0.1 Mq kpc^ Myr“^ and at the same 
time they do not keep a constant sign of Lx for the du¬ 
ration of the last 4 Gyr of the simulation. 

Eig.fTOl shows fractions of the three types of orbits as a 
function of radius in the final outputs of the simulations. 
The fraction of long-axis tubes correlates tightly with the 
amount of the mean rotation velocity of particles dis¬ 
cussed in section ICT Galaxies with higher rotation have 
faster increase of the fraction of long-axis tubes with ra¬ 
dius, mainly at the expense of the low angular momenta 
orbit family. 


5.3. Orbital frequencies 

Eor the long-axis tubes we computed the average an¬ 
gular and radial frequency of the stellar motion in the 
plane perpendicular to the major axis. The calculations 
were done for the period of the last 4 Gyr of each simu¬ 
lation. The angular frequency, cj, was computed in each 
time step as Vtjryz^ where is the tangential part of 
{Vy + and ryz = We defined the ra¬ 

dial frequency, as 1/Tavr, where Tavr is the average 
time period between two subsequent maxima of 

The distributions of the stars in the ujn plane for all 
four simulations are shown in the left panels of Eig. [TTl 
Only stellar particles with orbits classified as long-axis 
tubes are included. The right panels show the distri¬ 
bution of the ratio n/uj for the same particles. Two 
branches are clearly visible in the plots. Eor all simu¬ 
lations the branches have a similar slope of around 

1.8 for the upper branch and 1.4 for the lower one. Two 
sub-families are known to exist among the family of long- 
axis tube orbits in triaxial systems: the inner and outer 
long-axis tubes. A natural expectation therefore is that 
the two branches correspond to these two sub-families. 

In order to check this interpretation we tried to iden¬ 
tify t he inner and o uter l ong-axis tubes using the criteria 
from iGaida et al.l (j2Q15f ) . Eor an orbit to be classified 
as an inner tube, the absolute value of the x-coordinate 
of the particle when it reaches the maximum value of 
ryz must be larger than the corresponding value when 
ryz has a minimum, that is, \x{ryz,mEi^)\ > |:r(r^^,niin)|. 
Eor the outer tubes, the inverse condition is fulfilled: 

We find that indeed the outer tubes are more likely to 
be found in the upper branch of the k{uj) plot while the 
inner tubes are typically situated in the lower one. By 
visual inspection of some orbits we see that this rule does 
not identify the tubes correctly in all cases as particle or¬ 
bits are not always regular: they are still being disturbed 
by perturbations related to the past merger or by the 
time-dependent potential due to figure rotation. Eig. [12] 
shows two examples of orbits from the upper branch of 
k{uj) in the top row and from the lower branch in the 
bottom row that were classified as the outer and inner 
long-axis tubes respectively. Clearly, the shape of the 
orbits is as expected for these two sub-families, i.e. the 
outer tubes are convex while the inner ones are concave. 

Comparing the distributions of the ratio k/uj in the 






Fig. 11.— Left column: the distribution of the stars on long-axis 
tube orbits in the ouk (i.e. angular vs. radial frequency) plane. 
Right column: the distribution of the ratio The rows from 

top to bottom show results for simulations IGOR, I90R, I120R and 
I90VY1. 
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Fig. 12.— Examples of stellar orbits belonging to the long-axis 
tube family. The upper row shows three projections of the outer 
long-axis tube and the lower one of the inner long-axis tube. The 
orbits occupy different branches of the relation. 

right panels of Fig. [11] we see significant differences 
among simulations. In particular, the shape of the his¬ 
togram for simulation I60R is significantly different than 
for the other ones. The difference is not only that 
there are more stars on long-axis tubes (the distribu¬ 
tion takes larger values) but the left peak of n/uo ~ 1.4 
is significantly more populated than the right peak with 
n/uj ^ 1.8. Since we have identified the left peak (lower 
branch) with the inner long-axis tubes this means that 
the merger remnant in simulation IGOR is populated more 
by these orbits than the outer long-axis tubes. The sit¬ 
uation is opposite for the remaining simulations. Since 
the inner tubes tend to produce more boxy shapes this 
explains why the galaxy formed in IGOR was much more 
boxy (see Fig. 0 than the other ones. 


6. MERGERS VERSUS TIDAL STIRRING 

As mentioned in the Introduction, an alternative, and 
far better studied, scenario for the formation of dSph 
galaxies in the Local Group involves their interaction 
with a bigger host galaxy like the Milky Way or M31. 
In this section we attempt a comparison between our 
merger simulations and those describing such tidal stir¬ 
ring scenario. For this purpose we use four simulations of 
d warfs orbiting arou nd Milky Way-sized host described 
in lLokas et al.l (|2Q15f ). The initial model of the dwarfs in 
that paper was the same as we used for our merger pro¬ 
genitors and the individual simulations differed in the 
angle between the angular momentum of the dwarf disk 
and its orbital angular momentum. The four simulations 
will be referred to here by labels indicating the scenario 
(TS for tidal stirring) and the inclination angle: TS-IO, 
TS-I90, TS-I180 and TS-I120. TS-IO and TS-I180 cor¬ 
respond to the exactly prograde and exactly retrograde 
orbits, respectively, while TS-I90 and TS-I270 indicate 
the two cases where the dwarf’s disk was perpendicular 
to the orbit. 

As already discussed by iLokas et al.l (|2Q14al) , although 
using only an exactly prograde orientation of the disk, 
the main and most important difference between the two 
scenarios seems to be the occurrence of prolate rotation 
which is present in dwarfs formed by mergers and absent 
in the tidal stirring simulations. We are not aware of 
any mechanism which could induce prolate rotation only 
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Fig. 13.— Evolution of the mean rotation velocity around the ma¬ 
jor axis from simulations of the tidally stirred galaxies and Dwarf 
I of the merger simulation I90R. (Here, the colors used for other 
merger simulations throughout the paper were temporarily loaned 
to illustrate the results of the simulations of tidal stirring.) 


via tidal interaction of a disky dwarf with a host galaxy. 
To support this statement, in Fig. [13] we plot the evo¬ 
lution of the mean rotation velocity around the major 
axis, Vx^ of the stellar component of the dwarfs in the 
tidal stirring simulations. The measurements invariably 
yield very low values of this quantity, of the order of 1 
km s“^ at most, not only for the prograde case but for 
the other inclinations as well. As expected, the signal is 
strongest for the perpendicular orientations of the disk, 
i.e. in simulations TS-190 and TS-1270. We also note 
that the amount of prolate rotation is not related in any 
way to whether and how fast a bar forms and evolves 
in these simulations. For comparison, we also include in 
this Figure the same quantity measured for Dwarf 1 of 
our merger simulation I90R from Fig. [T] where the evolu¬ 
tion of Vx for all the merger simulations is shown. In the 
tidal stirring simulations, the rotation around the minor 
axis, Vz^ is always the dominant component of the rota¬ 
tion even though it is gradually diminished and replaced 
by random motions as a result of the evo lution (see the 
top panel of Figure 2 in lLokas et al.|[2QT5l) . 

Fig. [14] shows the surface brightness and kinematics 
maps for tidally stirred dwarfs from simulations TS-IO 
and TS-270 after G.5 Gyr of evolution (i.e. at the fourth 
apocenter) viewed along the intermediate and minor axis 
{y and respectively). The maps are analogous to those 
shown in Fig. [6] for the merger remnants, we kept the 
same range of the color scale as well as distances between 
the isocontours for the corresponding quantities. Only 
the displayed area is smaller here because the tidally 
stirred dwarfs are diminished in size at this late stage of 
evolution due to tidal stripping. The maps of the mean 
velocity (middle columns) clearly show that the rotation, 
if present, is seen always along the major axis of the pro¬ 
jected galaxy so the stars rotate around the minor axis. 
The velocity dispersion of the tidally stirred dwarfs (right 
columns of Fig. [T4|) is generally lower and more regular 
than in our merger remnants where the dispersion main¬ 
tains a more complicated structure even many Gyr after 
the merger. In order to detect these irregularities one 
would need a good spatial resolution for the kinematics 
and we comment on this issue further in the next section. 

The surface density maps for the tidally stirred dwarfs 
in the left columns of Fig. [14] show the remnant bar in 
the center. In comparison with these, the merger rem¬ 
nants have overall more regular shapes with more boxy 
isocontours. However, the boxiness is suppressed for non- 
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Fig. 14. — Maps of the surface density (left column), the mean 
line-of-sight velocity (middle column) and the line-of-sight velocity 
dispersion (right column) for simulations TS-IO and TS-I270 of 
tidally stirred galaxies. In the first and third row the line of sight 
is along the intermediate axis of the galaxy, in the second and 
fourth row it is along the minor axis. 

radial mergers and larger initial angles between angular 
momenta of the progenitor discs. In addition, the bar 
induced during the tidal stirring can make the isophotes 
look boxy as well if only the brightest central part of the 
galaxy was available for observation, so this feature prob¬ 
ably cannot be used as an indicator of the merger origin. 
In the case of low surface brightness galaxies, such as 
And II, the boxy shape would be particularly challeng¬ 
ing to detect. 

In principle, one more feature that could discriminate 
between the tidally stirred and merged objects, as adver¬ 
tised by iTomozein et al.l (j2Q15[) , is th e density profile in 
the o utskirts. It has been shown (e.g. iLokas et ani2Q13L 
l2Q15f ) that tidally stirred galaxies develop a flattening 
of the density profile due to the transition between the 
stars which are bound to the dwarf and the tidal tails. 
However, the presence and detectability of this transi¬ 
tion depends on the number of issues. Eirst, the transi¬ 
tion in the form of the break radius in the stellar density 
profile is not seen in dwarfs in non-prograde o rbital con¬ 
figurations (see Eigure 6 in iLokas et al.l 1^015f ). Second, 
the tails will not be visible if they are aligned with the 
line of sight of the observer and this unfortunately is the 


case for most of the time for dwarfs on diff erent orbits 
(|Klimentowski et ~aI]l2QQ9l: ILokas et al.ll2Q13[ ). Third, the 
tidal tails are usually very faint and have been convinc¬ 
ingly detected only for the Sgr dwarf. Moreover, tails 
produced by the merger can, at some stages, mimic this 
transition. Eurthermore, one cannot rule out the case of 
the merger remnants being subsequently tidally stirred 
by the host galaxy. 

7. TOWARDS A COMPLETE MODEL OF AND II 
EVOLUTION 

In this section we compare the results of our merger 
simulations to the real data for And II and discuss the 
possible scenario for its formation. In order to reproduce 
the observational data, we require the simulated merger 
remnant to show: (1) the rotation along the minor axis 
at the level of 10kms“^ at the maximum, which should 
be comparable to the maximum value of the line-of-sight 
velocity dispersion, (2) no significant rotation along the 
major axis, (3) the same kinematics for stars originating 
from both progenitors, (4) the shape of rather low ellip- 
ticity. The observational data cover the area of about 2 
kpc in radius. 

The measured kinematics of And II (|Ho et, a], I l2nn 
shows significant variability, while the equivalent values 
taken from the final outputs of the simulations are rather 
smooth if they are derived using all stellar particles along 
the major/minor projected axis about 4 kpc long. In this 
case even the fluctuations of the dispersion are averaged 
out for the final outputs. Therefore, in order to repro¬ 
duce the observed kinematics we explored the kinemat¬ 
ics from all simulation outputs (after the merger) as well 
as different lines of sight along which these outputs are 
‘observed’. Best matches of the observed kinematics (in 
terms of the lowest for all simulations are shown in 
Eig. [m In this comparison we again exclude the non- 
radial simulation 190VY2 which failed to produce signif¬ 
icant amount of prolate rotation. The sizes of the bins 
for the simulated data correspond to those used in obser¬ 
vations. The best matches come from times 4.6-7.1Gyr 
since the beginning of the simulations. In all cases we also 
verified that for the chosen line of sight the ellipticity of 
the galaxy i s in the range of about 0.1-0.2, as required 
by the data (iMcConnachie fc Irwinll2QQ6l: iHo et al.l 12 01 21: 
[Salomon et al.ll^l5f). 

In general, earlier stages of the merger can repro¬ 
duce the fluctuations in the measured kinematics bet¬ 
ter. Shortly after the merger, the galaxy is unrelaxed 
with many shells and streams visible that formed as a 
result of the collision. An example of the merger rem¬ 
nant image corresponding to the output at 4.7 Gyr with 
the best-fitting kinematics for simulation I90R is shown 
in the left panel of Eig. [TTl 

However, we are able to reproduce the fluctuations in 
the kinematics even for the final outputs corresponding 
to times a few Gyr after the merger if we recall that the 
real data are based on a much lower number of stars. Eor 
this comparison, we pick a subsample of about 500 stellar 
particles which match the n umber of s t ars us ed to derive 
the kinematics of And H in iHo et all (j2012[ ) . Again we 
use the same binning and select the stellar particles so 
that their distribution roughly covers the area used in 
the observations. In all other aspects, the choice of the 
particles is random. We repeat the subsample selection 
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Fig. 15.— Comparison between the observed (black points; 
IHo et ani2012l i and simulated (color lines) kinematics of And II. 
The panels show, from top to bottom, the mean line-of-sight ve¬ 
locity and line-of-sight velocity dispersion along the major axis of 
the observed galaxy (x) and the mean line-of-sight velocity and 
line-of-sight dispersion along the minor projected axis (y). For 
simulated galaxies, the line of sight is inclined with respect to the 
minor axis of the stellar ellipsoid by [27,167] deg for simulation 
I90VYI (4.6 Gyr; green lines), [18,166] for I60R (5.65 Gyr; cyan), 
[9,32] for I90R (4.7Gyr; red), and [32,158] for II20R (7.1 Gyr; ma¬ 
genta) . 

several hundred times for different lines of sight and find 
the ones that match the data best. 

Examples of good matches of the simulated data to 
the observed kinematics are shown in Fig. [16] for cases 
with ellipticity around 0.1-0.2. This proves that the ob¬ 
served fluctuations could be still consistent with the re¬ 
laxed stage of the merger remnant. An image of such 
a relaxed merger remnant in a late stage is shown in 
the right panel of Fig. [17] for simulation I90R at 10 Gyr 
where it can be directly compared to the analogous image 
at 4.7 Gyr in the left panel. These outputs and lines of 



y [kpc] 


Fig. 16.— The same as Fig. I15l but for a subsample of about 500 
stellar particles from the final outputs (at 10 Gyr) of the simula¬ 
tions. The line of sight is inclined with respect to the minor axis of 
the stellar ellipsoid by [45,127] deg for I90VYI, [18,114] for I60R, 
[27,126] for I90R, and [36,137] for I120R. 

sight correspond to the kinematics shown with red lines 
in Fig. [15] (for 4.7 Gyr) and Fig. [16] (for 10 Gyr). 

Besides the measured shape and kinematics, addi¬ 
tional constraints on the possible scenario for the for¬ 
mation of And II come from its recently measured star 
formation history (IVargas et al.l[^Ql4 IWeisz et al.ll2Ql4 

iGallart et al.ll2Ql^ . The measurements point to the ex¬ 

istence of the intermediate-age stellar population dated 
5-8 Gyr ago. Although we do not model star formation 
processes here, given the comparisons discussed above 
and these additional constraints we may speculate on 
the possible general scenario for the formation of And II. 

And II most likely formed as a result of a merger be¬ 
tween two disky gas-rich dwarf galaxies of similar mass 
and disk properties. The disks must be inclined by 90±30 
deg in order to possess significant components of angu- 
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Fig. 17.— Surface density of the stellar component in simulation 
is inclined with respect to the minor axis of the stellar ellipsoid by 
at 10 Gyr. 

lar momentum along the merger axis. The two galaxies 
approached each other on a nearly radial orbit with a 
rather low relative velocity forming a bound pair. The 
existence of such pairs of subhalos has been confirmed in 
recent studies of simulations reproducing environments 
similar to the Local Group (iKlimentowski et al.l [2QT(1 

Kazantzidis et al 1 l2Qlll:[Deason et al.ll2Ql4IWetzel et al l 

2Q15[ ). Since the occurrence of such groups is more likely 
outside the virial radius of big galaxies, the merger prob¬ 
ably occurred before the galaxies were accreted by M31, 
a few Gyr ago. The merger induced the secondary burst 
of star formation, mainly in the center of the merger rem¬ 
nant, converting some of the gas into stars in a similar 
fashion as in the well-studied me r gers between norma l- 
size galaxies (e.g. iCox et aT]l2QQ6l: iHoffman et al.ll^oTol ). 
The new stars gave rise to the distinct stellar popula¬ 
tion with a more concentrated density profile that can 
be identified with the intermediate-age stellar popula¬ 
tion seen in the data. These younger stars should con- 
tribute about 30% o f the stellar mass in the central region 
(jGallart et al.l[2Ql^ . 

Further evolution of the merger remnant remains an 
open question. However, since we do not see any young 
stars in And II and no gas is detected in it at present, we 
may suppose that the star formation ceased a few Gyr 
ago, either quenched by internal processes such as super¬ 
nova feedback or by environmental effects such as ram 
pressure stripping, if the merger remnant was accreted 
by M31. The present distance of And II from M31 is 
about 180 kpc and their relative line-of-sight velocity is 
of the order of 100 km s“^ i.e. low enough for And II to 
be a satellite of M31, but no strong constrain t s are avail- 
able on its orbit. However, iShava fc Tnllvl (j2013f ) find 
a solution for the Local Group galaxy orbits (see their 
Figure 13) where And H has experienced a close passage 
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)0R at 4.7 Gyr (left panel) and 10 Gyr (right panel). The line of sight 
1,32] deg for the output at 4.7 Gyr and by [27,126] deg for the output 

near M31 which may have stripped its gas. Nevertheless, 
this interaction could not be too strong so that And H 
could preserve its prolate rotation and probably did not 
lead to any new starburst. Gases of accretion of merger 
remnants by bigger hosts have been also confirmed in the 
simulat ions of the Local Group and Milky Way envir on- 
ments (|Kazantzidis et aTll2Qlll: [Tomozeiu et al.|[2Q15f ). 

The scenario we sketched here can be developed into 
a full model of a gas-rich merger only by incorporating 
gas dynamics and star formation. The outcome of such 
simulations will depend on a number of parameters like 
the threshold for star formation and the amount of feed¬ 
back. These parameters can in principle be adjusted to 
reproduce the observed mass fraction of the newer stellar 
population. We expect that such a model will be able to 
preserve the prolate rotation of the old stellar population 
but it remains to be seen what will be the kinematics of 
the younger stars. A more complete scenario for the for¬ 
mation of And H will be described in detail by S. Fouquet 
et al. (in preparation). 

8. SUMMARY AND GONGLUSIONS 

In this work we studied the origin of prolate rotation in 
dSph galaxies forming via mergers of two disky dwarfs. 
Our work was motivated by the disc overy of such type o f 
rotation in the dSph galaxy And H. iLokas et al.l (j2Q14a[ ) 
demonstrated that this kind of rotation can indeed re¬ 
sult from a merger between two disky dwarf galaxies. 
They proposed a detailed, but rather fine-tuned scenario 
with very special initial conditions. Here we showed that 
these initial conditions can be varied to some extent and 
can still reproduce the main observational properties of 
Andll. 

Masses of the luminous parts of the dwarfs were cho¬ 
sen so that their combined luminosity is similar to 
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what is measured for And 11. The dark matter compo¬ 
nents were adjusted so that the kinematics of the rem¬ 
nants m atches the observations. In addition. iLokas et al.l 
(|2Q14af ) showed that the different density profiles of the 
two stellar populations could be reproduced by adopting 
different initial scale-lengths of the disks of the dwarfs. 
We did not address this issue in the present paper and for 
simplicity assumed that the two progenitors are identical 
in terms of structural parameters and masses. 

The evolution of the three radial mergers we considered 
is straightforward to interpret. In these simulations the 
initial conditions are symmetrical for both dwarfs. The 
initial angle between the angular momenta of the dwarf 
disks determines how much of the disk rotation is trans¬ 
formed into the prolate rotation of the merger product. 
With the lower value of the angle, corresponding to the 
more vertical orientation of the disks, we get more prolate 
rotation as well as a more boxy shape of the surface den¬ 
sity of the galaxy. The dispersion and, more importantly, 
the rotation in simulation II20R is not high enough to 
match the data while I90R and I60R both reproduce the 
line-of-sight kinematics of And II satisfactorily. 

From non-radial simulations (190VYI and 190VY2) we 
see that the perpendicular component of the relative ve¬ 
locity, corresponding to the orbital angular momentum, 
cannot be increased too much. Both simulations have 
the radial component of the relative velocity of 8kms“^. 
In simulation I90VYI, with the perpendicular velocity 
component of I kms“^, the merger remnant still roughly 
meets the observation constraints (although with less ro¬ 
tation and smaller dispersion than observed). The rem¬ 
nant of simulation 190VY2, with the perpendicular ve¬ 
locity component of 2kms“^, does not reproduce the 
data well: it has very small rotation and the kinematics 
of stars originating from different dwarfs is significantly 
different. We note that for non-radial mergers the initial 
conditions were set up so that the merger is prograde 
for one progenitor and retrograde for the other one. In 
combination with a larger orbital angular momentum, 
this resulted in a very different evolution in simulation 
190VY2 compared to 190VYI. 

All remnants showing significant prolate rotation have 
similar orbital structure, dominated by box orbits in the 
center of the galaxy and by the long-axis tubes in the 
outer parts. However, there are subtle differences de¬ 
pending on the initial conditions: the transition from the 
box orbits to the long-axis tubes takes place at smaller 
radii for remnants with stronger prolate rotation. In ad¬ 
dition, the remnant with the highest rotation and the 
most boxy shape (I60R) turns out to possess a notice¬ 
ably bigger contribution of inner versus outer long-axis 
tubes than any other remnant. 

We conclude that prolate rotation may be produced in 
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mergers with a variety of initial conditions. The inclina¬ 
tion of the disks with respect to the orbit must be sig¬ 
nificant in order to provide enough angular momentum 
in the remnants but can be varied by a few tens of de¬ 
grees, producing remnants of slightly different kinematics 
and shape. Some orbital angular momentum during the 
merger, up to a factor of a few of the intrinsic angular 
momentum of the merging disks, is also allowed and does 
not destroy the symmetry between the two components 
of the remnant. 

We compared the predicted properties of our merger 
remnants to those characteristic of galaxies evolving in 
the vicinity of bigger host and affected by tidal stirring. 
Although a number of features, such as the boxiness of 
the shape, irregularities in the kinematics and strong 
tidal extensions, could possibly help us to distinguish 
the dSph galaxies formed in the two scenarios, they are 
generally not accessible with present observational tech¬ 
niques. We propose that the only reliable feature that 
can act as a fingerprint of the past merger is the pres¬ 
ence of prolate rotation such as the one observed in And 
II. Such type of rotation can not be induced via tidal 
stirring at a significant level even for a variety of orbital 
configurations. 

We have also compared our merger remnants to the 
real data available for And II. We find that in order to 
reproduce the data the merger could not have happened 
recently, but rather a few Gyr ago. This statement is 
supported by the fact that the surface brightness distri¬ 
bution in And II is rather regular without obvious distor¬ 
tions and the need to be consistent with the recently mea¬ 
sured star formation history of the galaxy which points to 
the presence of an intermediate-age, rather than a young, 
stellar population. The subsequent evolution of And II 
should have been rather quiescent in order to preserve 
prolate rotation, although it may have interacted with 
MSI, which led to the stripping of its gas. 
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